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Vibrational dephasing of nitrogen molecule is known to show highly interesting anomalies near 
its gas-liquid critical point. Here we present theoretical and computational studies of the Raman 
linewidth of nitrogen along the critical isochore. The linewidth is found to have a lambda shaped 
temperature dependence near the critical point. As observed in experimental studies, the calculated 
lineshape becomes Gaussian as the critical temperature (T c ) is approached. Both the present 
simulation and a mode coupling theory (MCT) analysis show that the slow decay of the enhanced 
density fluctuations near the critical point (CP), probed at the sub-picosecond time scales by 
vibrational frequency modulation, along with an enhanced vibration-rotation coupling, are the 
main causes of the observed anomalies. 

PACS numbers: 64.70.Fx, 82.53.Uv, 61. 25. Em, 78.30. Cp 



As the gas-liquid critical point of a fluid is approached, 
several dynamical properties of the system, such as the 
thermal conductivity, the bulk viscosity, the linewidth 
of the central Rayleigh spectrum, show anomalous be- 
haviour. This behaviour is collectively known as critical 
slowing down 0,0 and is physically explained by exploit- 
ing the divergence of the correlation length and using the 
dynamic mode coupling theory (MCT) 0, 0] . In the re- 
cent years, several anomalies have also been observed in 
molecular relaxation processes near the critical point or 
in the supercritical fluid. One of them is the surprising 
augmentation of the solvent density around a solute[j|. 
Another interesting case is the vibrational phase relax- 
ation which shows more than one anomaly, yet to be 
understood. 

In their experiments, Clouter et aZ.Q found that the 
isotropic Raman lineshape of a simple fluid like N2 ex- 
hibits a remarkable additional non-rotational broaden- 
ing near the gas-liquid critical points {p C rit,Tcrit)- They 
measured the Raman spectra along the triple point to 
the critical point and the behaviour of the lineshape as 
the critical point is approached from above. Recently 
Musso et aZ.pJ measured the temperature dependence of 
the lineshape parameters (i.e shift, width and asymme- 
try) both along the coexistence and the critical isochore 
of liquid nitrogen and found that the temperature depen- 
dent linewidth (T) is A shaped. The lineshape was found 
to undergo a change from Lorentzian (away form T c ) to 
Gaussian (near T c ). 

The observation of such anomalies has, till now, de- 
fied a convincing explanation. Mukamel, Stern and Ro- 
nis [8| had earlier interpreted the rapid broadening of 
Raman lineshape as a manifestation of dynamical criti- 
cal phenomena. This interpretation raises the following 



questions. As the microscopic time correlation function 
involved decays in less than a picosecond (actually in 
about 200 fs), why and how arc the slow long wavelength 
density fluctuations important and relevant? The physics 
here is clearly different from the one involved in the crit- 
ical slowing down at the long wavelengths. 

In a recent investigation, vibrational phase relaxation 
of the fundamental and the overtones of the N-N stretch 
in pure nitrogen was simulated by MD simulations, and 
the mode coupling theory (MCT)0 was used to explain 
the simulation results. This study could reproduce the 
rapid rise in dephasing rate as the critical point is ap- 
proached along the co-existence line, although it didn't 
include the vibrational coordinate (q) dependence of the 
inter-atomic potential and ignored the cross-term be- 
tween the vibration-rotation coupling and force and their 
derivatives. Everitt and Skinner [lfjj studied the Raman 
lineshape of nitrogen in a systematic way by including of 
the bond length dependence of the dispersion and re- 
pulsive force parameters along the coexistence line of 
nitrogen. They also incorporated the cross-correlation 
terms which were neglected earlier. Their results for 
the lincshifts and lincwidths along the gas-liquid coex- 
istence of N2 are in good agreement with experimental 
results. These theoretical studies did not consider de- 
phasing along the critical isochore. In addition, a con- 
vincing picture of anomalies did not emerge. The results 
presented in this Letter provide first microscopic expla- 
nation of the above anomalies. 

The theories of the vibrational dephasing are all 
based on Kubo's stochastic theory of the lineshape[Tl|. 
extended to the study of vibrational dephasing by 
Oxtoby 12]. The isotropic Raman lineshape, I(oj), is the 
Fourier transform of the normal coordinate time correla- 
tion function, Cq (t) ^3 > 
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A cumulant expansion of Eq.(l) followed by truncation 
after second order gives the following well-known expres- 
sion of Cq(*)H3], 



< Q(t)Q(0) > = Re exp(iaj Q t + i < Aoj > t) 



x exp 



dt'(t - t') < Aw(t')Aw(Q) > 



(2) 



The frequency modulation time correlation function, 
C u (t) =< Auj(t)Auj(0) >, derives contributions from the 
atom-atom(AA), resonance(Rs), vibration-rotation(VR) 
coupling, and also the cross-terms We calculated the 
linewidth, the lineshape and the dephasing time of N2 for 
different thermodynamical state points of nitrogen, both 
along the coexistence line and the critical isochore, using 
Eqs. 1 and 2. 

The Hamiltonian of homonuclear diatomic molecules 
can be written as 



(3) 



where H„ is the vibrational Hamiltonian, T(q) is the total 
translational and rotational kinetic energy, U(g) is the 
inter-molecular potential energy, and q is the collection 
of vibration coordinates qi . The inter- molecular potential 
energy is sum of the following site-site (%) between two 
molecules i and j ^(| • The total intermolecular potential 
energy is taken to be 
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where n a jp = \rjp(qj) - n a (qi)\ and fjp(q-j) = ^(O) + 
fjpqj/2. rjp{qj) is the position of the nucleus of atom 
(3 in molecule j, and fj^ is the unit vector oriented from 
the center-of-mass of molecule j to the /3th atom. The 
vibrational coordinate dependence of Lennard- Jones pa- 
rameters are given by, e iaj[ 3 = y/e ia ejp = e(l+jq l +jq j + 

2^ 2 q iqj ) and a iaj p = = a (l + % + <%). We 
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use the linear expansion coefficients 7 and 8 as deter- 
mined by Everitt and Skinner 0, . 

Microcanonical (NVE) MD simulations [TEj were car- 
ried along the coexistence line and along the critical iso- 
chore using the leap-frog algorithm for different thermo- 
dynamical state points of nitrogen. A system of 256 di- 
atomic particles was enclosed in a cubic box, and periodic 
boundary conditions were used. In simulation, the sys- 
tem was allowed to equilibrat e for 100 00 time steps with 
At = 0.0002t, where r [= yj (mcr 2 /e), m being the mass 
of the molecule ] is found to be equal to 3.147 ps. The 
averages were calculated over 400 000 MD steps. The 
thermodynamic state of the system can be expressed in 
terms of the reduced density of p* — pa 3 and a reduced 
temperature of T* = fcgT/e. The density of the system 
has been expressed in terms of number of molecules per 



unit volume times <r 3 and the temperature is in units of 
e/fes- Where a and ks is the Lennard- Jones diameter of 
the molecule and is the Boltzmann constant, respectively. 
Limited number of simulations have been done with N = 
512 molecules. We found no significant difference in the 
linewidth at the larger system. 




FIG. 1: The simulated lambda (A) shaped linewidth (F) of ni- 
trogen along the coexistence line (solid square) and the critical 
isochore (solid circle). The experimental results for linewidth 
along coexistence line (open square) and the critical isochore 
(open circle) reported by Clouter and Kiefte and by Chesnoy 
in Ref. are also shown. The CP is indicated by an arrow 
on the abscissa. 

Figure 1 shows both the simulated temperature depen- 
dent dephasing rate of nitrogen and also the experimen- 
tal one. The solid squares are the linewidth for differ- 
ent state points along the coexistence curve of nitrogen 
and the solid circles are the linewidth along the criti- 
cal isochore. The interesting feature in the figure is the 
lambda-shaped linewidth when the values for two dif- 
ferent regions of nitrogen are presented together. This 
figure is similar to the one observed experimentally (see 
figure 4 of Ref.0)- It is interesting to note the sharp rise 
in the dephasing rate as the CP is approached. There 
are noticeable difference in the high temperature region 
along the critical isochore. 

To understand the origin of this critical behaviour, 
we carefully analyzed each of the six terms |l6l| (three 
auto correlations and three cross-terms between density, 
vibration-rotation coupling and resonance) , which are re- 
sponsible for the modulation of the vibrational frequency 
for fundamental transition. Two terms are found to 
dominate near the CP and these are the density and 
the vibration-rotation coupling. The temperature depen- 
dence of these two terms are shown in figure 2 where the 
integrand of Eq.2, X y -(i) = J* dt' (t-t')C z J (?), is plotted 
against time for seven state points along the critical iso- 
chore. Where i, j represents the density (den), vibration- 
rotation coupling(VR) and resonance(Rs) terms respec- 
tively. Note the sharp rise in the value of the integrand 
as the critical temperature is approached, and the fall 
when it is crossed. We have found that both these con- 
tributions at the CP are distinct compared to the other 
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state points. Thus, the rise and fall of dephasing rate 
arises partly from the rise and fall in the density and the 
vibration- rotation terms. 
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FIG. 2: The time dependence of the relative contributions 
of the (a) density-density, Xden(t), and (b) VR coupling, 
Xvn{t), terms at different temperatures along the critical iso- 
chore. 

A crossover from Lorentzian-like to Gaussian lineshape 
can happen when the usually large separation in the time 
scales of decay of C w (t) and Cq (t) 0, H3 ceases to exist 
and the two time correlation functions begin to overlap. 
Indeed, the computed lineshape becomes Gaussian near 
the from CP. Note that the frequency modulation time 
correlation function decays fully in about 200 femtosec- 
ond. 

What could be the reason for this dramatic crossover 
behaviour? The increase in density fluctuations 0, 0] 
near the CP increases the magnitude of the mean square 
frequency fluctuation < Aujf >. This leads itself to an 
increase in the rate of decay of the normal coordinate 
time correlation function, (Cq(£)). In addition, the fre- 
quency modulation time correlation function slows down 
somewhat. As the decay of these two time correlation 
functions become comparable, the lineshape goes over to 
the Gaussian form. This a novel effect and different from 
that commonly encountered. 

If one considers only the number density as the rele- 
vant slow variable in dephasing, then mode coupling the- 
ory analysis (MCT) gives the following expression for the 
density dependence frequency modulation time correla- 
tion function [l8||. 

(A<j p (0)Au p (t)) = / k 2 dkF s {k,t)C 2 {k)F(k,t)^ 



where C{k) is the Fourier Transform of two particle di- 
rect correlation function. Near CP, the main contribu- 
tion is derived from the long wavelength (that is small 
k) region, where self-intermediate scattering function, 
F 8 (k,t) = e~ Dsk * and intermediate scattering function, 
F(k, t) = S(k)e- DTk2t , D s is the self-diffusion coefficient, 
where S(k) is the static structure factor, and Dt is the 
thermal diffusivity. Thus, (Aw p (0)Au p (t)) ~ S(k 
0) e -(D,+D T )k 2 t_ Near cp g( k _> o) becomes very large 
( as compressibility diverges at T = T c ), leading to- 
wards a Gaussian behaviour for lineshape. Dt also un- 
dergoes a slowdown near T c . This complex dependence 
may lead to a Levy distribution from time dependence 
of (Q(0)Q(t)) as discussed earlier by Mukamel et al.Q 
However, a limitation of the above analysis is the absence 
of the VR term which contributes significantly and may 
mask some of the critical effects. At high temperature, 
the latter dominates over the density term. A complete 
Lorentzian behaviour is predicted only in the low tem- 
perature liquid phase. Interestingly, the predicted diver- 
gence of < Aw 2 (0) > very close to T c enhances the rate 
of dephasing and this shifts the decay of < Q(t)Q(0) > 
to short times, giving rise to the Gaussian behaviour. 
We shall return to this point later. We have used MCT 
to demonstrate that the large enhancement of Vibration- 
Rotation coupling near the gas-liquid critical point arises 
from the non-Gaussian behavior of density fluctuation 
and this enters through a non-zero value of the triplet 
direct correlation function. 

To further explore the origin of these anomalous criti- 
cal temperature effects, we have investigated for the pres- 
ence of dynamical heterogeneities in the fluid at three 
temperatures near the CP, by calculating the well-known 
non-Gaussian parameter a(t) defined as|l9(, a(t) = 

(|) ^ 2 nr^i — 1 where < Ar(t) 2 > is the mean squared 

displacement and < Ar(i) 4 > the mean quartic displace- 
ment of the center of mass of nitrogen molecule. It can 
only approach zero (and hence Gaussian behaviour) for 
times exceeding the time scale required for individual 
particles to sample their complete kinetic environments. 
As can be seen from figure 3, the function a(t) is large 
near CP at times 0.5 - 5 ps, indicating the presence of 
long lived heterogeneities near T c . 

The contribution of the density term reflects the com- 
bined effects of collisions, repulsive and attractive inter- 
actions on the friction on the bond. This term naturally 
grows as the CP is approached. The rise in the contribu- 
tion of the vibration-rotation term has a different origin. 
This increase is due to the slowing down of the rota- 
tional time correlation function near the CP [2(3. As we 
approach the CP, the number of molecules having slow 
decay of rotational time correlation function increases. 
While one can describe the average effect of density fluc- 
tuation through the divergence of the static structure 
factor S(k) at small wavenumbers (as discussed above), 
this approach does not capture the full scenario. This is 
because the decay of the frequency modulation time cor- 
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FIG. 3: The non-Gaussian parameter a(t) is plotted against 
time (t) at three different temperatures along the critical iso- 
chore. 

relation function (C u (i)) occurs in the femtosecond time 
scale. At that time scale, the density fluctuation is nearly 
static. Thus, vibrational dcphasing provides snap shots 
of the large density fluctuations present near the critical 
temperature. Even more interesting is the origin of the 



ultrafast decay of C^(t)-this is partly due to the cancel- 
lation which arises from the cross-terms of VR coupling 
with the density and the resonance terms. These terms 
have sign opposite to the pure terms but have compa- 
rable magnitude, leading to further enhancement in the 
rate of decay of C w (i). The reason for such large negative 
cross-correlation can be understood in terms of the inho- 
mogeneity. Note that critical anomaly is nearly absent 
in CF4, CH4 and CC^Hl). While cubic symmetry in the 
the former two precludes VR coupling, C0 2 may be too 
heavy for VR to be important. 

It is indeed surprising that our simulations could cap- 
ture many of the novel features observed in experiments, 
including the lambda shaped temperature dependence of 
the dephasing rate and the cross-over from the Lorentzian 
to the Gaussian form. Large density fluctuations near T c 
shifts the dynamics probed to sub-picosecond times. This 
combined with VR coupling, give rise to the observed 
anomalies. 
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